Non-exponential relaxation and hierarchically constrained dynamics in a protein 
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A scaling analysis within a model of hierarchically constrained dynamics is shown to reproduce the 
main features of non-exponential relaxation observed in kinetic studies of carbonmonoxymyoglobin. 
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There have been various works in which similarities 
between the dynamics of proteins and the structure and 
dynamics of glasses and spin glasses have been discussed 
111- LJ- Li Ljl • Although the energy landscapes in the two 
systems may be quite different, non-exponential relax- 
ation is observed in both. In glasses, relaxation often 
follows the stretched exponential form characteristic of 
the "Kohlrausch law:" 



*(t) =$ exp[-(t/r) /3 ], Q<0<1, 



(1) 



where r is a temperature-dependent characteristic time 
which becomes unmeasurably long as the glass transition 
temperature is approached. It is often experimentally 
observed to follow, over 10 orders of magnitude, a Vogel- 
Fulcher law 0, r oc exp[^/(T - T )] 

Now, any reasonable relaxation function $(t) can be fit 
by assuming some distribution w(t) of relaxation times 
among additive contributions to the relaxing quantity, 
thus 



*(<) 



dTw(r) exp(— t/r). 



(2) 



This extends the idea of conventional Debye relaxation 
with a single relaxation time to a situation in which there 
is a distribution of degrees of freedom each contributing 
independently to $(i) with its own relaxation time - thus 
parallel relaxation. 

A different point of view was proposed by Palmer, 
Stein, Abrahams, and Anderson |fg. They pointed out 
that the conventional parallel picture, while simple, is 
often microscopically arbitrary and that a more physi- 
cal view is that the path to equilibrium is governed by 
many sequential correlated steps - thus a series inter- 
pretation in which there are strong correlations between 
different degrees of freedom. These authors proposed 
a microscopically motivated model of hierarchically con- 
strained dynamics (HCD) which leads to the Kohlrausch 
law (and a maximal relaxation time of the Vogel-Fulcher 
form). That result was cited by Shibata, Kurata, and 
Kushida Q to argue that HCD holds in their observation 
of stretched exponential relaxation in an experiment on 
conformational dynamics in Zn-substituted myoglobin. 
Of course it is possible, and sometimes likely, that both 
parallel and sequential processes exist in the same sys- 
tem. 
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FIG. 1: Experimental relaxation in MbCO Q 



In a discussion of anomalous relaxation in proteins, 
Iben et al |3j presented data on carbonmonoxymyoglobin 
(MbCO). This is myoglobin in which CO is bound to the 
central iron atom of the heme group. In this system, par- 
allel relaxation processes also occur; they dominate the 
dynamics of the rebinding of the ligand after photodisso- 
ciation Q • Here we focus on the pressure release relax- 
ation experiments of Iben et al. 3] Infrared absorption 
spectra of the stretch bands of the CO were taken under 
various conditions. After pressure release, the center fre- 
quency of the Aq band initially shifts rapidly upward by 
0.4 cm -1 and then relaxes slowly toward its low-pressure 
equilibrium value, 1.2 cm -1 higher. This behavior at var- 
ious temperatures is shown in Fig. 3 of Ref. 1 and it is 
rcplotted here (without the error bars) as Fig. 1. It seen 
that the relaxation is close to power law over more than 
three decades of time, thus much slower than exponen- 
tial, or even stretched exponential. 

It is of interest to ask whether there is a model of hi- 
erarchically constrained dynamics as proposed in Ref. 
which can account for the main features of Fig. 1. These 
are: a) At each temperature there is a region of power-law 
relaxation which crosses over at shorter times to some- 
thing much slower, b) This crossover time increases as 
temperature decreases, c) The power law is the same at 
each temperature, but there is an increasing offset as the 
temperature is lowered. 
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In what follows, it will be shown that one of the mod- 
els mentioned in Ref . |6| in fact gives behavior identical to 
what was observed in MbCO. For HCD, one recognizes 
that equilibrium distributions in configuration space are 
not relevant since the free energy barriers which deter- 
mine relaxation are continuously changing in time as dif- 
ferent degrees of freedom relax at different rates. Fur- 
thermore, in a strongly correlated system one expects 
that with any choice of coordinates, interactions will re- 
main in the form of constraints and that these will be of 
importance over a range of time scales. The nature of 
constraints is that some degrees of freedom cannot relax 
until their motion is made possible by the relaxation of 
other degrees of freedom. These restrictions occur over 
a wide hierarchy of coordinates, from fast ones to slow 
ones. A complete discussion of this HCD approach is 
given in Ref. |y. 

To implement this picture, Palmer et al. |(| set up a 
hierarchy of levels 71 = 0,1,2,3,.... The degrees of free- 
dom in level n are represented by N n Ising pseudospins 
(two-level centers) each of which has two possible states. 
This was adapted from earlier work of Stein 0. Con- 
straints enter via the requirement that each "spin" in 
level 7i + 1 is free to change its state only if a condition 
on some number /i„ of spins in level n is satisfied. Now, 
the fin spins have 2^™ states. Let the required condition 
be that just one of these possible states is realized. If the 
average relaxation time in level n is r n , then on average, 
it will take a time 2^"r n for a spin in level n+1 to change 
its state. Therefore 

n-l 

r n = t exp(/Ji In 2) = r exp U n , (3) 
o 

where 

n-l 


The relaxation function is given as a sum of the correla- 
tion functions of all the degrees of freedom Sf. 

N 

#(t) = (l/i\0£<Si(0)Si(t)>. (5) 
«=i 

In a correlated system, the dynamics of the Si are not in- 
dependent, so each of the correlation functions in Eq. (5) 
depends on the behavior of the other S"s. As described 
above, the HCD scheme of Palmer et al. [|| incorporates 
correlations. The Si are arranged in a hierarchy of levels 
7i with each level having its characteristic relaxation time 
r n given by Eq. (3). So Eq. (5) may be rewritten as a 
sum over the different levels from to oo: 

oo 

$(t) = ^™ n exp(-i/0, (6) 

n=0 



where w n = N n /N is the fraction of the total number of 
degrees of freedom which are in level n. 

For a given model, w n and /i n must be specified. As 
remarked in Ref. |g, the simple choices fi n — /Uo/hi2, a 
constant and w n = w n -i/X give power-law relaxation. 
Here, this situation is examined more fully. 

The sum in Eq. (6) is rewritten as an integral using 
the above choices for w and ji: 

$(t)=w [ dn\- n exp{(-t/T )e- n ^}, (7) 
Jo 

where from normalization at t = 0, wq = In A. This 
integral is evaluated exactly in terms of the incomplete 
gamma function 

7(1/, u) = / dxx v - Y e- x . (8) 
Jo 

The result is 

Ht)=^(jYlMro), (9) 

where v = In X/fio . For large values of its second ar- 
gument, j(v,u) approaches the complete gamma func- 
tion T(u). Thus, at large times $(t) cx t~ v , as observed 
3]. For small u, j{v,u) « u u /v - u v+l /{v + 1) + . . ., 
so that for short times, $(t) crosses over to a slower 
1 — const(i/ro) dependence. Temperature dependence 
is introduced in the model through the temperature de- 
pendence of the fundamental relaxation time tq(T); its 
inverse corresponds to the rate constant k r (T) introduced 
in Ref. |3- Thus the behavior of Eq. (9) is similar to the 
form 

*(t) = [l + k r (T)t]- v (10) 

which was used in Fig. 3 of Ref. to fit the data. 

The appearance of the experimental points at different 
temperatures, in particular that the data are parallel (the 
same power law for all T) at long times, suggests that a 
scaling function could describe the experimental results. 
A general form is 

<f(t,T)^T a G[t/T (T)}, (11) 

Therefore, by rescaling the time by a parameter To(T) 
for each temperature the data for $>(t,T)/T a would all 
fall on a single curve. The fact that the power law of 
the long time behavior is independent of T implies that 
the exponent a = 0. If this rescaling is carried out, the 
temperature dependence of the characteristic rescaling 
time tq(T) may be determined. This rescaling for the 
data of Iben et al Q is shown in Fig. 2a. 

If in addition one has a theoretical expression which 
has a scaling form, such as Eq. (9), then by collapsing the 
data onto the theoretical curve, one obtains the numerical 
values of tq(T). This is carried out for the hierarchical 
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FIG. 2: (a) Rescaled data of Iben et al. |3j (b) Comparison 
of theory and rescaled data 

model of Eq. (9) in Fig. 2b. The value of v in Eq. (9) is 
adjusted to agree with the common large t slope as seen 
in Fig. 2a. The black curve in Fig. 2b is Eq. (9) with v 
set equal to 0.28. The values of To for each of the four 
temperatures are plotted in Fig. 3, where the solid curve 
is a fit using the expression 

l/T- (T) = fc exp[-(To/T) 2 ] (12) 

The result is k = 1.5 x 10 16 5s" 1 , T = 1105 K.[ljj 
With only four points, a number of different functional 
forms might seem equally good. In particular, an Arrhc- 
nius form works almost as well , bu t has an unacceptably 
large pre-exponential factor, [llj The form chosen has 
been suggested by several authors 0, ; it can describe 
diffusion in a random potential of a form which mimics 
the potential surface of a protein. 

As discussed in Ref. 0, the observed relaxation process 
most likely involves substates having different conforma- 
tional molecular structures as well as different angles of 
the CO ligand with respect to the heme normal. With 
this in mind, contact may be made between the parame- 
ters of the above model and the actual experiment. Re- 
call that no measures the number of "spins" in a given 
level which must arrive at a certain configuration before 
a typical spin in the next level becomes unfrozen and can 
relax. The physical picture of constrained movement of 
atoms which underlies the model and its application to 
MbCO then implies that this number should be around 
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FIG. 3: Fundamental relaxation rate 



3 to 5 or so. The number A measures the geometric re- 
duction of the fraction of spins which belong to level n 
as n increases through slower and slower levels. $(£) has 
only a logarithmic dependence on A so it is reasonable to 
take In A to be of order unity. This gives v — In X/fiQ i$ 1- 
The fit value is v = 0.28, i.e. slightly less than four spins, 
on average, combine to unfreeze the degrees of freedom 
at the next level. This seems quite reasonable. 

Eq. (10), taken from Eq. (1) of Ref. |3j is consistent 
with the present Eq. (9) in that they agree in the lim- 
its of small and large i/ro- Therefore the fit by Iben et 
al. 0] using Eq. (10) is also satisfactory. However, there 
is no physical motivation for that form whereas in the 
present work the observed non-exponential relaxation in 
MbCO is derived from a microscopic model which incor- 
porates strong-correlation constraints on the relaxation 
of the molecule in a definite way. 

What conclusions can be drawn from the present anal- 
ysis? The scenario that the model describes is one in 
which the primary relaxation event represents the rate 
1/tq(T) at which a typical enthalpy barrier is over- 
come. Since no further T-dependence is introduced, the 
conclusion is that all subsequent relaxation events are 
"slaved" to the primary one and that they represent 
entropic conformational changes of the molecule. 

A Los Alamos group has independently been analyzing 
the relaxation processes in MbCO. [ill fl2l | Remarkably, 
they have reached conclusions which are consistent with 
the present hierarchical model for the relaxation of the 
CO vibration frequency. Namely, they argue that the 
temperature dependence of the relaxation is governed by 
an activation enthalpy for which the solvent is respon- 
sible. This determines the rate of the fastest relaxation 
process - 1/tq(T) in the hierarchical model. Subsequent 
relaxations involve degrees of freedom of the protein and 
the hydration shell and are governed by entropy barri- 
ers; they have the same temperature dependence as the 
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solvent fluctuation rate and are slower. This descrip- 
tion is precisely the same as that of the hierarchical 
model presented here. The physically motivated scenario 
of the successful hierarchical approach lends support to 
the identification of the physical relaxation processes de- 
scribed in Ref. 12. 

The analysis presented here can be generalized to more 
complicated situations. For example, the hierarchical 
rules could be modified to include simultaneous paral- 
lel relaxation ( "unslaved" ) processes as in the ligand re- 
binding referred to earlier, internal enthalpic barriers, 
reverse constraints, and intra-level correlations. While 
other forms than Eq. (9) might fit the MbCO data, within 
the hierarchical scenario the fact that the long time be- 
havior is a scalable power law practically forces the simple 
rules which were used to obtain Eq. (9). The success of 
this approach suggests that a similar picture and analysis 
can be relevant for other dynamical processes in biolog- 
ical molecules. If so, insight can be obtained about the 
physical processes which determine the relaxation phe- 
nomena. 
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